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Abstract: We investigate the large- iV behavior of the topological susceptibility x 
four-dimensional SU(A^) gauge theories at finite temperature, and in particular across 
the finite-temperature transition at T^. For this purpose, we consider the lattice 
formulation of the SU(A^) gauge theories and perform Monte Carlo simulations for 
= 4, 6. The results indicate that x has a nonvanishing large- iV limit for T < Tc, 
as at T = 0, and that the topological properties remain substantially unchanged 
in the low-temperature phase. On the other hand, above the deconfinement phase 
transition, x shows a large suppression. The comparison between the data for = 4 
and = 6 hints at a vanishing large- limit for T > Tc. 
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The pattern of chiral symmetry breaking for QCD with Nj hght flavors at zero 
temperature is well understood. The symmetry group of the classical lagrangian, 
U(iVj)L ® \J{Nf)ji, is broken both by the anomaly and spontaneously The sponta- 
neous breaking of the axial subgroup SU(iVj)^ yields N"^ — 1 Goldstone bosons, while 
the anomalous breaking of the U(1)a symmetry explains the heavier fiavor-singlet 
state observed in the hadronic spectrum. (See e.g. Ref. [|T|.) 

In the large- limit 0, where N is the number of colors, the mass of the singlet, 
m^/, is related to the topological susceptibility of the pure gauge theory, x, through 
the well-known Witten-Veneziano (WV) formula [§, |^ : 

F>2, = 2NfX. (1) 

It is particularly interesting to study the topological susceptibility as is varied, 
since, at fixed number of flavors, 1/N can be identified with the explicit symmetry- 
breaking parameter for the U(l)^ symmetry. Hence QCD is expected to recover 
the full U(A^/)l ® U(A^/)r chiral symmetry as oo. It is then possible to show 

that, under very general assumptions, this chiral symmetry is spontaneously broken 
down to the vector subgroup \J{Nf)v, yielding Nj massless Goldstone bosons p. In 
this limit, if the topological susceptibility of the pure gauge theory does not vanish, 
the rj' acquires a mass m^, ~ i.e. it becomes a Goldstone boson whose mass 

squared vanishes linearly in the symmetry-breaking parameter 1/iV as the anomaly 
is suppressed, see e.g. [^. At zero temperature, numerical evidence from lattice 
simulations in favor of a no n- vanishing large- A^ limit of with 1/A^^ power-law 
corrections, has only been obtained recently [0, |]. 

As the temperature is increased, the SU(A^)a chiral symmetry is restored at a 
critical temperature T^. The nature of this phase transition, besides its own theo- 
retical interest, determines the dynamics of the transition from hadronic matter to 
a quark-gluon plasma, which is expected to take place, e.g. in heavy-ion collisions. 
In this respect, the effective breaking of the U(1)a symmetry at finite temperature, 
and in particular around the transition at T^, is of particular interest. Indeed, in 
the case of two light flavors, the transition may be continuous (for massless quarks) 
and in the 0(4) universality class only for a sufficiently large breaking of the U(1)a 
symmetry around 0, [T^. At finite temperature the anomaly, considered as an 



equation between operators, remains unchanged [|TT|, but its physical effect might be 
drastically different because the quantities that enter the WV formula have their own 
temperature dependence. In order to verify that the WV mechanism is still at work 
at finite temperature and in the low-temperature phase, we study the behaviour of 
the topological susceptibility x the critical temperature is approached from be- 
low, verifying that x has a nonvanishing large- A^ limit, similarly to what happens 
at T = 0. As a consequence, Eq. (P is expected to hold up to Tc. The transition 
between the high- and low-temperature regimes is not fully understood. At high 
temperatures, where instanton calculus is reliable [IT^, a rather different scenario 
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emerges. Concerning the topological properties, an interesting hypothesis has been 
put forward in Ref. |T^: at large N, configurations with non-trivial topological charge 



are exponentially suppressed in the high-temperature phase, i.e. behave as e~^, so 
that the topological susceptibility gets rapidly suppressed in the large- limit. In 
order to shed some light on this issue, we also present results from simulations above 
T 

The behavior of the topological susceptibility across the transition has already 
been investigated for N = 2 and = 3 in a number of works, see e.g. Refs. |]14|, 



15, 16, 17, Qq, M. The main focus of this work is on the large- behavior. For this 



purpose we report numerical results for SU(A^) gauge theories with = 4, 6. Some 



results for A^ > 3 have recently been reported in Ref. |T9|. 

The behavior of the topological susceptibihty at the finite-temperature decon- 
finement phase transition is studied in detail for SU(A^) gauge theories with A^ = 4, 6, 
exploiting their lattice formulation given by the action 

S = -N^ Yl Tr [U^{x)U^{x + + i^Mix) + h.c] , (2) 

X,fl>U 

where U^{x) G SU(A^) are link variables. In order to study the large- A^ limit, it is 
convenient to replace the more familiar coupling /5 by 7 = /3/2A^^, which is the inverse 
of the 't Hooft coupling A = g'^N . In order to study the theory at finite temperature, 
we perform Monte Carlo simulation on asymmetric lattices. The gauge configurations 
are generated using a mixture of microcanonical and heat-bath updating algorithms 
(see Ref. pO| , pl| for details). In our simulations we consider different time extensions 



Lt = 6,8 and constant aspect ratio Lg/Lt = 4. These values of Lt and Ls should 
be sufficiently large to obtain results with small scaling and finite-size corrections. 
The physical temperature is given as a function of the lattice spacing and of the 
lattice time extension, T = l/a{'y)Lt. Previous investigations have already shown 
the existence of a finite-temperature phase transition, which is first order for A^ > 3, 



see e.g. Ref. ||2^. For each value of Lt, 7 is tuned so as to explore an interval 



around the critical temperature Tc] the corresponding critical value of the coupling 



is denoted by 7c. In this work the critical couplings obtained in Ref. are used, i.e. 

= 6) = 0.33717(2) and jc{Lt = 8) = 0.34640(7) for A^ = 4, and -f^{Lt = 6) = 
0.34508(5) for A^ = 6. When zero-temperature quantities are needed, we perform 
the corresponding computation on symmetric lattices at the required values of the 
coupling. A summary of our runs is presented in Tables |l| and ^, for A^ = 4 and 
A^ = 6 respectively. 

The range of 7 values that can actually be explored with the Wilson action is 
rather limited. On the one hand, 7 has to be sufficiently large so that the system is 
actually in the weak-coupling region, i.e. beyond the first-order bulk phase transition 
at 7 = 0.339 in the case A^ = 6, and beyond the crossover region characterized by a 
peak of the specific heat at /5 = 0.325 for A^ = 4; see Ref. PO, I2TI] for a more detailed 
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Table 1: Finite-temperature data for the SU(4) gauge theory. ot=o and xt=o are re- 
spectively the string tension and the topological susceptibility at T = obtained using 
symmetric lattices. Data marked by an asterisk are subject to an uncontrolled system- 
atic error due to the fact that no clear plateau was observed in the cooling procedure to 
determine x- 



discussion of this point. On the other hand, as 7 is increased, two types of difficulties 
arise. First, as the lattice spacing becomes smaller, larger lattices are necessary to 
avoid finite-size effects. In practice, we always try to use values of 7 such that the 
spatial extent satisfies Ls\^ > 3, where a is the zero-temperature string tension, 
as suggested by previous investigations H, |20|, ^ The second obstacle is the 
increase of the autocorrelation time of the topological modes as the continuum limit 
is approached |23[. Such a severe form of critical slowing down puts a stringent 
limit on the upper value of 7 that can be efficiently simulated, at least with the 
currently available algorithms. Finally, since the transition is first order and rather 
strong in the case N = 6, some attention must be paid to the dependence of the 
Monte Carlo results on the starting configurations. In particular, close to the finite- 
temperature phase transition, we use hot or cold starting configurations according 
to the side of the transition we were investigating, to avoid hysteresis effects. 

We use the following formulas and definitions for the rescaled and reduced tem- 
peratures: 

Tr{Lt,i)^^= i^ , (3) 



,ir ,_ T,(L„7)-T,(L„7c(Lt)) hMLt)) 

t{Lt, 7) = — — — = W 1, (4) 

Tr{Lt,-fc{Lt)) Y ct(7) 

where, as before, a is always computed on symmetric lattices, i.e. at T = 0. In order 
to determine the reduced temperature according to Eq. (H), we used the following 
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Table 2: Finite-temperature data for the SU(6) gauge theory. crT=o and xt=o are re- 
spectively the string tension and the topological susceptibility at T = obtained using 
symmetric lattices. 



values of the string tension at 7c: y/a = 0.270(2) for 7 = 0.33717 and = 4 
(obtained on a 12^ x 24 lattice), ^/a = 0.203(1) for 7 = 0.3464 and = 4 (obtained 
on a 16^ X 32 lattice), and y/a = 0.276(2) for 7 = 0.34508 and = 6 (on a 12^ x 24 
lattice). 

The topological charge Q is estimated using the cooling technique described in 
Ref. 0]. We compute the corresponding susceptibility x = (Q^)/^ the scaling 
ratio: 

R{LuV = —, r- 5 

x(oo,7) 

The results are reported in Tables |] and |^ respectively for = 4 and A^ = 6. 

Before discussing the results for the topological susceptibility, let us assess the 
possible sources of systematic errors in the lattice computation. It is well known that 
topological structures may disappear during the cooling procedure. This problem is 
particularly severe at finite temperature, especially for small Lf, where one may 
not observe clear plateaux with respect to the cooling steps, see e.g. Ref. [^. At 
zero temperature, a direct comparison with a fermionic estimator of the topological 



charge shows good agreement ||25|, supporting the idea that the coohng method 
is fairly stable in this case. However, the situation at finite temperature is more 
difficult to control and can eventually generate a systematic error. In this respect, 
we expect the cooling method to perform better with increasing A^. In order to 
minimize the possible bias due to the cooling technique, only lattices with Lt > 6 
have been considered in this study. In Fig. ^we show measurements of the topological 
susceptibility during cooling for A^ = 4, = 6, 8 and N = 6, Lt = 6, and values 
of the couplings corresponding to the high-temperature phase, and in particular 
t ~ 0.05. Plateaux during cooling are clearly observed (from 4 to 20 cooling steps) in 
the simulations for A^ = 6, for A^ = 4 with Lt = 8, and only in the low-temperature 
phase in the case A^ = 4 with Lt = 6. This allowed us to unambiguously determine 
the topological susceptibility using the cooling method also at finite temperature. On 
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Figure 1: Topological susceptibility along the cooling process. The values of 7 correspond 
to reduced temperatures t ~ 0.05 in all cases. 



the other hand, for = 4 in the high-temperature regime, the time-direction size 
Lt = 6 was not sufficient to provide a clear plateau, as shown in Fig. |l| for the data 
corresponding to t ~ 0.05. Measurements based on the cooling method become rather 
questionable in these cases, because they do not provide an unambigous estimator of 
X- In Table |I| the corresponding results are marked by an asterisk: they were obtained 
after 10-12 cooling steps, but they are subject to an uncontrolled systematic error, 
which probably leads to an underestimate of x, and they should therefore be taken 
into account only with due care. For small N and Lt other estimators for Q should 
be used. 

The data for the scaling ratio R are displayed in Fig. |^. One can immediately 
remark that its behavior is drastically different in the low- and high-temperature 
phases. In the low-temperature phase, all data for = 4, Lf = 6, 8 and = 6, 
Lf = 6 appear to lie on the same curve, showing that scaling corrections are small 
and also that the large- limit is quickly approached. The ratio R remains constant 
and compatible with 1.0. Only close to Tc, i.e. for T > 0.97 Tc, does this ratio 
appear to decrease. These results show that in the confined phase the topological 
properties remain substantially unchanged up to T^. On the other hand, above the 
deconfinement phase transition, we observe a large suppression of x- The comparison 
between the A^ = 4 and A^ = 6 data shows that the ratio R decreases much faster 
for A^ = 6, hinting at a vanishing large- A^ limit of R for T > Tc. 

Numerical results supporting the same picture have also been reported recently in 



Ref. [T^]. The numerical evidence of the topological suppression across the transition 
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Figure 2: The ratio i? as a function of the reduced temperature t. 



was inferred from simulations at T^, by monitoring the correlation of the topologi- 
cal charge with the Polyakov line, whose value is used to infer the actual phase of 
the configurations generated along the given Monte Carlo run. Therefore a more 
quantitative comparison with our results is not straightforward. A comparison with 
the results presented in Refs. |T^, |TB[ suggests that the suppression of topological 
fluctuations is faster in SU(4) than it is in SU(3). 

A nonvanishing topological susceptibility x implies a nontrivial dependence on 
the 6 term that appears in the euclidean Lagrangian as 



4 



(6) 



where q{x) is the topological charge density. Indeed x is the second derivative of 
the free-energy density (ground-state energy) F{9) with respect to 6' at 6^ = 0. More 
generally, expanding the free-energy density around 6 = 0, one may write 

1 



F{0) = -xO' (l + M' + M' + ...) 



(7) 



The parameters of the expansion of F{6) are related to the moments of the probability 
distribution P{Q) of the topological charge Q in the large- volume limit. While x is 
determined from the second moment {Q"^), the coefficients 62* are related to higher 
moments of P{Q), for example 
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The large-volume limit of the probability distribution P{Q) is Gaussian only if b2i = 
0, i.e. when F{9) = without higher-order corrections. A nontrivial expansion 
around 6 = 0, such as Eq. (|^, reflects deviations from a simple Gaussian behavior 
of P{Q), whose size is controlled by the coefficients • ^ The coefficient 62 has been 
estimated in Ref. for the SU(A^) gauge theories with = 3, 4, 6 at T = 0, finding 
very small values, i.e. 62 = -0.023(7) for = 3, 62 = -0.013(7) for A^ = 4, and 
62 = -0.01(2) for A^ = 6, supporting the conjecture ^ 62 = 0(1/ N^). Thus 
for A^ > 3 the simple Gaussian form F{9) ^ is expected to provide a good 

approximation of the dependence on 6 for a relatively large range of values of 6. In 
order to investigate this issue at finite temperature, we have also computed 62 in 
our finite-temperature Monte Carlo simulations. In the low-temperature phase the 
estimates of 62 turn out to be compatible with those at T = 0, suggesting that F{9) 
remains substantially unchanged up to T = T^, with very small corrections to the 
Gaussian behavior. On the other hand, in the high-temperature phase the absolute 
value of 62 turns out to be significantly larger, for example at t ~ 0.05 we found 
62 — —0.05 for A^ = 4 and 62 — —0.08 for A^ = 6, indicating larger deviations from 
the Gaussian behavior, although they are still moderately small. 

In conclusion, the results presented in this paper suggest that the physical prop- 
erties determined by the fiuctuations of the topological charge, such as the Witten- 
Veneziano relation (|l]), remain substantially unchanged in the low-temperature con- 
fined phase. In particular, the topological susceptibility of the pure gauge theory 
has a nonvanishing large- A^ limit in the low-temperature phase, as at T = 0. On the 
other hand, in the high-temperature phase there is a sharp change of regime where 
the topological susceptibility is largely suppressed. Such suppression becomes larger 
with increasing A^, suggesting that the topological charge vanishes above the critical 
temperature. Monte Carlo results seem to support the scenario presented in |TB[: at 
large A^ the topological properties in the high-temperature phase, for T > Tc, are 
essentially determined by instantons from very high temperature down to T^; the 
exponential suppression of instantons induces the rapid decrease of the topological 
activity observed in the large- A^ limit. 



^An apparently contradictory result has been reported in Refs. |2^, for the expected large- 

volume probability distribution P(Q), i.e. P{Q) = {2Tr{Q'^)y^^'^e ^(q^ [1 + 0(l/y)], which was 
obtained starting from a generic expansion of F{6) at 9 = 0, such as Eq. (|^), and evaluating the 
large-volume behavior of P{Q) = J d9 ex-p[—iQ6 — VF{d)] using a saddle point approximation. 
A Gaussian behavior in the large-volume limit [p7| would contradict the assumption of a generic 
expansion of F(6), and, in particular, a nonzero value of 62, which implies a nonzero fourth moment 
of the large- volume P{Q)- The point is that the contributions considered as 0{1/V) corrections 
to the Gaussian behavior cannot be neglected in order to reproduce the correct physically relevant 
large- volume limit of the distribution's moments. This can be checked by computing the corrections 
to the saddle-point approximation. 
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